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ABSTRACT 

Observations of the Lyman-a forest and of high-redshift galaxies at z ~ 6- 10 imply that there were just 
enough photons to maintain the universe in an ionized state at z ~ 5-6, indicating a "photon-starved" end 
to reionization. The ionizing emissivity must have been larger at earlier times in order to yield the extended 
reionization history implied by the electron scattering optical depth constraint from WMAP. Here we address 
the possibility that a faint population of galaxies with host halo masses of ~ 10 8 ~ 9 M Q dominated the ionizing 
photon budget at redshifts z k, 9, due to their much higher escape fractions. Such faint, early galaxies, would not 
have formed in ionized regions due to suppression by heating from the UV background, and would therefore not 
contribute to the ionizing background at z < 6, after reionization is complete. Our model matches: (1) the low 
escape fractions observed for high-redshift galaxies, (2) the WMAP constraint of r es ~ 0.09, (3) the low values 
for the UVB at z < 6, and (4) the observed star formation rate density inferred from Lyman-break galaxies. 
A top heavy IMF from Pop III stars is not required in this scenario. We compare our model to recent ones in 
the literature that were forced to introduce an escape fraction that increases strongly towards high redshift, and 
show that a similar evolution occurs naturally if low mass galaxies possess high escape fractions. 
Subject headings: cosmology: theory — dark ages, reionization, first stars — intergalactic medium 



1. INTRODUCTION 

A fundamental question in the theory of reionization is: 
what were the ionizing sources? While it was realized 
several decades ago that the observed quasar population 
was insufficient ( Sha piro & Girouxll 19871) . the actual sources 
have eluded detection. The main constrai nts on reioniza- 
tion m odels are the Gunn-Peterson trough dGunn & Peterson! 
119651) in the spectra of high redshi f t quasars at a red- 
shift of z ~ 6 (e.g., iFan et all 12001 iWillottl I2007T) . and 
the electron scattering optical depth, r es , with the latest 
results indicating that the universe w as substantially ion- 
ized by z ~ 10 (Koma tsu et al.l 12011). Models adhering 
to the se constraints predict ex t ended reionization histories 
(e.g. lHaiman & Holder! [2003] ICenl 12001 IWvithe & Loebl 



Fan et al 



2003; 



Ricotti & Ostriker 
20061: 



2004; Choudhurv & Ferrara 2006; 
Haiman&Brvan 2006: llliev et al.l 12007: 



Kuhlen & Fauch er-Giguerd 120121; lHaardt & Madaul 120121; 
Ahn etal 1120 121) . 

The ionizing emissivity o bserved at z ^4-6 is t oo low to ac- 
count for the observed r es . iBolton & Haehneil ((2007) found 
that extrpolating the ionizing emissivity inferred from the 
Lyman-a forest at z ^4-6 to higher redshifts would underpro- 
duce r es . Likewise, observations of "Lyman-break galaxies" 
(LBGs) at z > 6 indicate a rapidly declining sta r formation 
rate in the brightest galaxies toward high redshift dOesch et alj 
120121) . co nsistent with evo lution of the dark-matter halo mass 
function ( Trent feTaT] 120 1 Oh . The ionizing emissivity from 
these galaxies can reionize the universe by z ~ 6, under fairly 
generous assumptions regarding uncertainties such as the ion- 
izing escape fraction, the intergalactic medium temperature, 
and the abundance of galaxies fainter than the current de- 
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tection limit dTrenti et alj|2010b iBouwens et alj|201 ll) . Thus, 
both observations of the high-redshift galaxy population and 
of the Lyman-a forest indicate a photon-starved end to reion- 
ization. 

Reconciling the ionizing emissivity observed at z ~ 4 - 6 
with the t rs value has be e n the subject of recent work. 
iKuhlen & Faucher-Giguerd d2012l) presented a model that 
matches both a photon-starved end to reionization and the 
WMAP result. In their best-fitting model, the escape frac- 
tion, / esc , varies with redshift, while the relation between halo 
mass and ionizing luminosity is extrapolated from abundance 
matching at lower redshift. The "min imal reionization model" 
adopted by lHaardt & Madaul d2012l - HM12 hereafter) also 
assumes / e sc increases with redshift. iFontanot et al.l d2012l) 
find that either an increasing / esc to wards high redshift and/or 
towards fainter galaxies is required. IShull et aTl d2012l) argue 
instead for an end of reionization of z ~ 7 and a r es at the 
low-end of the WMAP allowed range. 

Here we investigate the scenario in which low mass galax- 
ies have a much higher / esc than more massive ones, but do 
not survive past reionizion. We separate halos into two cat- 
egories based on their mass. Low-mass halos are subject to 
complete suppression by a photoionizing background and are 
called "photosuppressible", while more massive halos are not 
sensitive to an ionizing background. Under these circum- 
stances, photosuppr ession effects naturally extend the reion- 
ization process (e.g.. lHaiman & Holderl2003llWvithe & Loebl 
l2^03llChoudhury & Ferrarall2006l: llliev et alJl2007h . We fur- 
ther assume that low-mass halos have similar star formation 
efficiencies to their high-mass counterparts, but have much 
higher escape fractions, of order unity. Before proceeding, we 
first discuss the two main physical motivations for our model. 

First, theoretical models and observations indicate the ion- 
izing photon / esc should vary with halo mass. Star for- 
mation in the lowest-mass galaxies was likely bursty, tak- 
ing place in a highly irregular gas distributions more re- 
sembling massive star forming regions than rotationally sup- 
ported disks. As illustrated by the numerical simulations pre- 
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sented by IWise & Cenl ( 120091) . stars should be able to carve 
out channels through which radiation can escape. For higher 
mass halos, the escape fraction should only reach a few per 
cent, because most of the ionizing radiati on is produced deep 
within the galaxy's disk and absorbed dGnedin et al.1 120081: 
iRazoumov & Sommer-Larsenll2010l) . In support of this pic- 
ture, observations at z ~ 3 indicate that Lyman Break galaxies 
possess low ionizing esc ape fractions (< 10%: IShaplev et all 
l2006t ISiana et ail 120101) while fainter gal axies show escape 
fractions that gro w comparable to unity (llwata et al.1 12009b 
iNestor et alll201 ll) . 

Second, photoionization heating has a negative feedbac k 
effect on star forma tion in dwarf galaxies (Efstathiou 1992). 
Shapi roet al.l (Q994) investigated the linear growth of gas per- 
turbations, finding a filtering scale that lags beh ind the in- 
stantaneous Jeans scale, and Gnedin & Hui ( 1998) presented 
convenient analytical formulae for the filtering scale. Sub- 
sequent coupled radiation-hydrodynamical simulations con- 
firmed the negative feedback effect of photoionization heating 
dThout & Weinberg 1996; lGnedinll2000l: iDiikstra et al.ll2004t 
Okamoto et al. 2008). The latest simulations bv lFinlator et all 
d201 ll) indicate a suppression scale of a few times 10 9 M Q . 
Photoionization heating will have a significant negative effect 
on star formation in the faintest, most abundant, early galax- 
ies. 

In this paper we attempt to find the simplest possible reion- 
ization model that satisfies current constraints. An important 
simplifying assumption we have made is the same halo mass 
threshold for the escape fraction and photosuppression - in 
our fiducial model, halos with masses below 2 x 10 9 Mq have 
both high escape fractions and are subject to photosuppres- 
sion. A common mass scale is physically plausible, however, 
given that it is the same mechanism (photoionization feed- 
back) which is likely to be responsible both for higher escape 
fractions and suppression in low-mass halos. 

We present our model in §2, results in §3, and a dis- 
cussion in §4. All distances and densities are in comov- 
ing units, and we adopt cosmologica l parameters consis - 
tent with the latest WMAP constraints dKomatsu et al .11201 ll) . 
(Cl m , A , n b , h, cr 8 , n s ) = (0.275 , 0.725 , 0.046, 0.7, 0.82, 0.97). 

2. MODEL 

We obtain the hydrogen ioniz ed fraction, x(z), using 
the following relation ship (e.g., I Shapiro & Girouxl 11987b 
lHaardt&Madaull2012l) : 



«0 (free) ' 

where n 7 is the comoving ionizing photon emissivity, tig = 
« H + « He is the comoving number density of hydrogen and 
helium nuclei, and the recombination time is (f rec ) _1 = 
»oQflCj(z)( l +z) 3 . We use the clumping factor c? = Cioo ob- 
tained by Pawliketal. (2009) for their z le j on ~ 10.5 case, 
C/ = e -° 28z+3 - 59 +l forz > 10andc/ = 3.2forz< 10, consistent 
with the transmissivity and mean free path reported in the lit- 
eratur e dSongaila & Cowidl2010b iKuhlen & Faucher-~G iguere 
[2012b . Helium is singly ionized along with hydrogen, with 
He II reionization occurring later, at z = 3. 

The ionizing emissivity is related to the star formation rate 
density (SFRD) via 

«7 =/esc,l/7P*,l+/esc.2/7/0*,2 (2) 

where f esc j corresponds to type i halos, / 7 is the number of 



ionizing photons produced per stellar mass formed, and p\j is 
the SFRD of type i halos. Type 1 (photosuppressible) halos 
can only host sources of ionizing radiation if they are in neu- 
tral regions, while more massive type 2 halos can host sources 
even after reionization. 

We use a Salpeter IMF, with M rain = 0. 1M Q , M max = 100M©. 
At solar metallicity wi th this IMF, / 7 ~ 3 x lO^M" 1 (e.g., 
lHaardt & Madaul [20121 and references therein), while for a 
metal poor population with Z = 0.02Z o , / 7 ~5x IO^Mq 1 , 
and for the Z = case / 7 ~ 7 x IO^Mq 1 (ISchaeredl2002l) . We 
assume that the ionizing photons are emitted instantaneously, 
which is an excellent approximation given the lifetime of the 
massive stars that produce the bulk of the ionizing radiation. 

The SFRD is given by 

p\ = p*,i+p*,2 = Po^p [(l-*)e*,i/i + e*,2/2] , (3) 

where e* ; is the fraction of baryons converted to stars in halos 
in the mass range i with collapsed fraction /, . The ( 1 —x) factor 
in the first term accounts for the suppression of star formation 
in ionized regions for photosuppressible halos. We will use 
Mi and M2 to indicate the lower mass limits of our two halo 
mass ranges. Equations ([T), (|2), and OJ are solved simultane- 
ously to obtain the ionized fraction, ionizing emissivity, and 
star formation rate density as a function of redshift. 
The electron scattering optical depth from reionization is 

3H n b ca T [°° x{z){\+zf{\-Y+N He {z)Y /A) 

T es = / ; dz, (4) 

where Y ~ 0.24 is the helium abundance and Afae(z) is the 
number of times helium is ionized - we use A^e = 2 for z < 3 
andA^He = 1 forz > 3. 

3. RESULTS 

Shown in Figure Q] are our results for our three models. 
The first two, "suppression" and "no suppression", were ob- 
tained using the model in §2. In both cases, Mi = 10 8 M Q 
(roughly the mass at which atomic cooling is efficient), M2 = 
2 x 1O 9 M (the mass at above which star formation can pro- 
ceed in ionized regions), / esc ,2 = 0.05, and 1 = e*.2 = 0.03. 
In the no suppression model, we removed the 1 — x factor on 
the right hand side of equation (01. We restricted e*^ < e* 2, 
since internal feedback effects are stronger at lower masses, 
and e*,2 = 0.03 was chosen so as to match the overall effi- 
ciency found of halos with M > 2 x 1O 9 M in the lTrenti et alj 
(2010) ICLF model (see Figure 1 and discussion below). We 
chose /esc,2 = 0.05 to be consistent with the relatively low es- 
cape fractions inferred for observed LBGs. We varied the re- 
maining parameter, / esc .i, in both cases to obtain r es = 0.086, 
so that /esc.i =0.8 and / esc ,i = 0.45 in the cases with and with- 
out suppression, respectively. 

The last model shown is our implementation of the "mini- 
mal reionization" model of HM12, with a constant value of 
/ 7 . In their implementation, the metallicity is assumed to 
follow Z(z) = Zq10~°' , which corresponds to a value of 
/ 7 ~ 5 x 10 60 Mq at z ~ 11, decreasing at lower redshifts. 
Since such an evolution is degenerate with a slight change 
/esc(z), we do not expect this to be important. We used 
their fit for p*(z), which was extrapolated beyond z ~ 7 from 
the observed data (overproducing p\(z =10) according to 
lOesch et aT]l2012l) . and an evolving escape fraction, f esc (z) = 



Reionization by the First Galaxies 



3 



i i i | i i i | i i i | i i i I i 

Suppression 



1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 
No Suppression 




8 10 12 
z 



FIG. 1 . — top: Ionized fraction x(z) and neutral fraction 1 —x(z)- The dashed line shows the mean escape fraction middle: Comoving ionizing photon emissivity 
(solid), as well as the contribut ion from halos in each of our two m ass ranges (dashed). Solid points with errorbars show measurements from the Lyman- 
alpha forest at z =5-6, as compiled by Kuhlen & Faucher-Giguere (2012, see their Table 2 and references therein). Dotted lines indicate the emissivity required to 
keep the universe ionized, ri 7 = «o/('rec), and are only plotted for redshifts such that x > 0.9. bottom: Comoving SFRD, p*(z), for each of the models (solid), as 
well as the contribution from halos in each of our two mass ranges (dashed). The triangular point with errorbars is the value obtained from GRB observations at 
z ~ 9 by Robertson & Ellis (2012). Square points indicate a fit to the observed SFRD for halos with M > 2 X 10 9 Mq, and its extrapolation to z > 8 iTr enti et alj 
l20TOh . 



1.8 x 10" 4 (1 +z) 3 - 4 . For the HM12 model, we also obtained 
r es = 0.086. In all three cases we used / 7 = 4.5 x lO^'M" 1 . 

3.1. With suppression 

As seen in panel (a) of Figure 1, the ionized fraction x 
passes 50% at z ~ 10 and 90% at ~ 7. Also shown in panel 
(a) is the mean emission-weighted escape fraction, defined by 



(/esc)(z) = 



fyP*(z) 



(5) 



The escape fraction rises from ~ 0.05 at z < 6, to ~ 0.8 at 
z = 1 1 . This evolution is entirely due to the evolution in the 
underlying source population - as the universe becomes more 
neutral, more ionizing photons come from photosuppressible 
halos, which have much higher escape fractions. 

The ionizing emissivity is sh own in panel (d), as well 
as con straints from table 2 of iKuhlen & Faucher-G iguere 
d2012l) . derived by combining measurements of the mean free 



path (e.g., Songaila & Cowie 2010) with the transm issivity of 
the Lyman-a forest (e.g., iBolton & Haehnem2007l) . It reaches 
a peak of n 7 ~ 10 51 s _1 Mpc~ 3 at z ~ 11. The upturn at z ~ 5 
is due to the rise of higher mass sources, in particular from 
halos above 2 x 10 9 M Q . The main sources of reionization 
in this model do not contribute substantially to the emissiv- 
ity at z < 6, consistent with both an early reionization that 
satisfies the WMAP optical depth, and a photon-starved end 
to reionization. The star formation rate density rises from 
p\ ~ lO" 3 M0yr _1 Mpc" 3 at z = 16, plateaus at ~ 0.01 over 
Z ~ 8— 12, increasing at z < 1 2. (panel g) . Show n also are 
values from the ICLF model of ITrenti et all ([20 10) integrated 
down to a minimum halo mass of 2 x 10 9 M Q - our model 
matches the predicted evolution of the SFRD of Lyman-break 
galaxies. By z ~ 13, star formation occurring in suppressible 
halos accounts for over 90 per cent of the total SFRD. 



4 



ALVAREZ, FINLATOR & TRENTI 



3.2. No suppression 

As seen in panel (b), overlap occurs earlier in the model 
without suppression, at z ~ 8.5 (as compared to z ~ 5.5 in 
the suppression model). The post-reionization ionizing emis- 
sivity is higher without suppression, n 7 > 10 51 s _1 Mpc~ 3 for 
z > 5 (panel e). This violates inferences from the Lyman- 
a forest at z ~ 5 - 6, which indicate values at those red- 
shifts of ~ 2-4 x 10 50 s~ 1 Mpc~ 3 (data points in Figure [p 
iBolton & Haehnemi2007t IKuhlen & Faucher-Giguerell2012l) . 
and would be difficult to reconcile wi th the appearance of the 
Gunn -Peterson trough at z 6 (e.g., iFan et al.1l2002t IWillottl 
l200l . We were unable to find a model with e* \ < e* j an d 
/esc,2 < 0.05 that simultaneously satisfies the constraints and 
does not include suppression. 

3.3. Minimal reionization model ofHM12 

In Figure 1, we show results from our implementation of 
the "minimal reionization model" of HM12. Shown in panel 
(c) is the escape fraction, / esc oc (1 + z) 3,4 . Since n 7 oc / eS cP*, 
the drop in the SFRD towards higher redshift in panel (i) is 
roughly canceled by the increasing / esc , so that « 7 is nearly 
flat (panel f). Notice that scenarios where the ionizing emis- 
sivity is flat struggle to reproduce the upper limit to the emis- 
sivity at z = 6. By contrast, the suppression model naturally 
accommodates this limit with an increasing contribution from 
fainter sources at higher redshifts (panel d). However, the up- 
per limit at z = 6 (« 7 < 2.6 x 10 50 s _1 Mpc~ 3 ) could be higher, 
as it is based on combin ing inferred photoionization rates 
( IBolton & Haehn elt 2007) with highly-uncertain e stimates of 
the mean free path at z = 6 (Songaila & Cowie 2010). A lower 
mean free path than assumed by [Kuhlen & Faucher-Giguere 
(120121) . which was ba sed upon a fit which lies above the 
z = 5.7 data point from Songaila & Cowie (2010), would re- 
sult in a higher upper limit. 

Both the HM12 model and our suppression model exhibit 
a similarly strong evolution in the mean escape fraction, with 
/esc of order unity for z > 10, and / esc ~ 0.05 at z ~ 5. In 
the suppression model, this is a natural outcome of the as- 
sumption that the low-mass halos have high escape fractions 
and dominate the ionizing photon budget early (/ eS c,i = 0.8), 
but by the end of reionization only the more massive halos 
host stars, and the escape fraction is about a factor of 15 
lower (/esc,2 = 0.05). Suppression can therefore be viewed as a 
physical explanation for the strongly evolving escape fraction 
needed to satisfy the observational constraints. 

Comparing the SFRDs assumed by the HM12 and the sup- 
pression model reveals that they agree to within a factor of 2 
for all z =5-13. In detail, the SFRD in the HM12 model uses 
an extrapolation to z > 8 of a fitting function constrained at 
z<8: 



6.9 xl0" 3 + 0.14 (z/2 ' 2) 



1.5 



l+(z/2.7) 41 



M yr _1 Mpc" 3 . (6) 



The suppression model predicts that the mass range dominat- 
ing star formation transitions smoothly from photosuppress- 
ible halos to more massive ones around z =8-10. In fact, com- 
paring the dashed curves in panel (g) with the SFRD in panel 
(i), the HM12 model depends on either the increasing role of 
sources that cannot be detected in the post-reionization Uni- 
verse (including photosuppressible halos), or a strongly in- 
creasing star formation efficiency in LBGs. 

4. DISCUSSION AND OBSERVATIONAL IMPLICATIONS 



In this letter we presented a reionization model that can si- 
multaneously satisfy all the current observational constraints 
on hydrogen reionization (Lyman-a forest at z ~ 5 - 6, SFRD 
from Lyman-break galaxy surveys at z ~ 6 — 10, low es- 
cape fractions at z ^ 6, and WMAP r cs ). Our model 
predicts that the majority of the ionizing photons at z > 
9 are produced in faint galaxies hosted in photosuppress- 
ible dark matter ha los with M ~ 10 8 5 undetectable 
in current surveys dTrenti et al.ll2010t iMunoz & Loebll2~01 It 
IKuhlen & Faucher-G iguere 2012). Once reionization is com- 
pleted, the star formation in these halos is suppressed by ra- 
diative feedback, so that these sources do not contribute to the 
ionizing emissivity at lower redshift. Therefore, we predict 
that the sources of reionization are primarily galaxies with 
L < 0.06L*(z = 3), too faint to be detected by even the deepest 
HST observations. 

A key prediction is that the SFRD is higher than that 
inferred from Lyman-break galaxy surveys, by a factor of 
~ 2 at z = 9, increasing to a factor of ~ 10 at z = 12. In- 
deed, alternate tracers of star formation indicate that this is 
likely. Estimates of from the GRB rate measure p\ ~ 
c" 3 at z - 9 dRobertson & Ellisll2012t 



Q.025; ^! 8 , M^yr'M ) 



yr 'M pi 

see also Kistler et all 12009), one ordei of magnitude higher 
than the LBG surveys estimate. In addition, the non-detection 
of GRB host galaxies at z > 5 in ultradeep HST images pro- 
vides an independent confirmation that the majority of star 
formation at those redshifts is h appening in galaxi es fainter 
than current observational limits ( Trent ietalJl20 1 2|) . Finally, 
the latest determination of the galaxy luminosity function at 
z ~ 8 has a very steep faint-end slope (a = -1.98 ± 0.2) at 
Mab ~ -17.7 , implying a logarithmically diver gent contribu- 
tion t o the SFRD from fainter, unseen galaxies (Bradley et al. 
120121) . All these observations are naturally explained in the 
framework presented here based on photosuppressible halos 
providing most of the reionizing photons. 

It is useful to compare the SFRD and escape fractions we 
obtain to that of the HM12 model, which is based on the 
SFRD in galaxies with L > 0.06L* at z = 2-7. In our frame- 
work the sources of reionization at z ~ 10 will remain unseen 
in near-future galaxy surveys at z ~ 10 - our photosuppress- 
ible halos correspond to -15 £s Mab ^ — 10 over z — 6-12. 
If the SFRD in the HM12 model comes from galaxies of the 
same UV luminosities as at lower redshifts, their extrapola- 
tion to z = 10 is likely too high and marginally inconsistent 
with the strong decline in the LBG popu lation between z = 8 
and z = 10 observed by Oesc h et all d2012l see their Fig. 8). 

Even assuming a SFRD due to bright galaxies that is not de- 
clining as quickly as suggested by the LBG observations, it is 
still necessary to assume an escape fraction that is a strongly 
increasing function of redshift in order to satisfy the WMAP 
r es constraint. Such an evolution is naturally explained in our 
model, where the average escape fraction increases to high 
redshift simply because there are more low-mass halos form- 
ing in the unheated regions. 

The high escape fractions for faint galaxies favored in 
our model should h ave a unique signature in the cosmic in- 
frared background dFernandez et al.l 120121) . This is because 
the Lyman-a emission per escaping ionizing photon is in- 
versely proportional to the escape fraction. A larger escape 
fraction for a fixed ionization history will result in a lower 
overall flux and fluctuation amplitude. 

There are several important caveats. First, we have as- 
sumed that the star formation efficiency in photosuppress- 
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ibl e halos is t h e sam e as more massive ones. Simulations 
by IWise et aTl d2012l) that include the effects of radiation 
pressure, supernova feedback, and metal-line cooling, find 
M, ~2x 10" 2 M Q yr _1 for a 2 x 10 8 M Q halo at z = 8, imply- 
ing that ~ 0.03 of the gas accreted after reaching the atomic 
cooling limit is converted to stars, consistent with our ch oice 
of e* ,i = 0.03. The simulations of iFinlator et al.l (1201 ll) . on 
the other hand, exhibit a much lower star formation rate at 
a similar mass of a few times 10 8 M Q . Although theoretical 
predictions have not converged, e*,i ~ e*,2 seems plausible. 

In addition, we have assumed that the suppression oper- 
ates instantaneously within H II regions. Although this is 
an oversimplification, it illustrates the basic trends resulting 
from high escape fractions and suppression in low-mass halos. 
Detailed numerical simulations that self-consistently include 
both external and internal radiative feedback during reioniza- 
tion will be required to address these issues more accurately. 

Finally, we note that the best-fitting value for the electron 



scattering optical depth from WMAR r es = 0.088, has 1-tr er- 
ror bars of of ~ 0.015, indicating that the true value could be 
as low as r es ~ 0.06, at 2-a. In our model with suppression, 
t ~ 0.06 can be achieved with a lower escape fraction for low 
mass halos, / esc .i ~ 0.2. This is still however a factor of four 
times higher than that for the larger mass halos, / esc = 0.05. 
Upcoming measurements with the Planck satellite will be es- 
sential in tightening constraints on r es and therefore on the 
nature of the first galaxies. 
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